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We present a study of the pentaquark system in quenched lattice QCD using diquark-diquark 
and kaon-nucleon local and smeared interpolating fields. We examine the volume dependence of the 
spectral weights of local correlators on lattices of size 16® x 32, 24® x 32 and 32® x 64 at /3 = 6.0. We 
find that a reliable evaluation of the volume dependence of the spectral weights requires accurate 
determination of the correlators at large time separations. Our main result from the spectral weight 
analysis in the pentaquark system is that within our variational basis and statistics we can not 
exclude a pentaquark resonance. However our data also do not allow a clear identification of a 
pentaquark state since only the spectral weights of the lowest state can be determined to sufficient 
accuracy to test for volume dependence. In the negative parity channel the mass extracted for this 
state is very close to the KN threshold whereas in the positive parity channel is about 60% above. 
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I. INTRODUCTION 

Experimental searches for a pentaquark state near the 
KN threshold, reported in recent experiments HQ , are 
under way world wide. The accumulation of evidence 
from low energy experiments for the existence of this 
state combined with the negative results obtained in 
high energy experiments 0 pose deep questions as to its 
nature and its production mechanism. Further doubts as 
to the existence of the O'*" have been raised since the re¬ 
sent report of the CLAS collaboration for lack of evidence 
for a resonance state from a dedicated high statistics pro¬ 
ton experiment |^. The 0+ was predicted theoretically 
in the chiral soliton model as an exotic baryon state 
with an unusually narrow width. The possible existence 
of such a state has raised interesting questions on what 
its structure should be in order to account for its nar¬ 
row width. A number of phenomenological models have 
been put forward to explain its stability such as special 
flux tube configurations Hiilli and diquark forma¬ 
tion [m. 

Several studies in lattice QCD have looked for a pen¬ 
taquark state in order to determine its mass and par¬ 
ity but no consensus has been reached yet with some 
groups finding a bound state with mass close to the ex¬ 
perimental value I 12L 1 131 jl4j. HI in and others the KN 
scattering state [ 13 , H mi’ One main difference be¬ 
tween these groups has been the interpolating field used 
to create this state [20) . Since its structure is unknown, 
optimizing the interpolating field that one uses is a dif¬ 
ficult task. However, for reasonable interpolating fields, 
the results should not be dependent on the interpolat¬ 
ing field. In this work we use two of the most popular 
choices: an interpolating field motivated by the diquark- 
diquark picture and an interoolating field motivated 
by the kaon-nucleon structure Q]. We note that diquark 
formation was shown to be important in the context of 
the static pentaquark potential where arrangements of 


quarks that allow diquark formation are found to be en¬ 
ergetically favored giving rise to a potential that is pro¬ 
portional to the minimal flux tube connecting the five 
quarks mmiiii Using these interpolating fields as 
a basis we construct a 2 x 2 mass correlation matrix 
and determine the optimal linear combination that yields 
maximum overlap with the ground state of the system. 
Another way of enhancing the ground state by suppress¬ 
ing excited state contributions is smearing. In addition 
to using local interpolating fields, we also use smeared 
fields by applying gauge invariant smearing to the quark 
fields used. We find that the value we obtain for the 
mass is independent of which of these various interpo¬ 
lating fields we use. Since both interpolating fields have 
an overlap with the KN scattering state we expect, that 
being the lowest state, it will determine the large time 
dependence of both correlators. If a pentaquark bound 
state exists about 100 MeV higher than the KN thresh¬ 
old we expect it to dominate the time dependence of the 
correlator up to time separations of about ~ 10 GeV“®. 
Therefore the main difficulty is to identify unambiguously 
the resonance in an intermediate time range before the 
nearby KN scattering state dominates. If our interpolat¬ 
ing fields have good overlap with the KN scattering state 
and the 0’*', then diagonalization of the mass correlation 
matrix should give as lowest eigenvalues these two states. 
This is our motivation for choosing the KN interpolating 
field, which is expected to have large overlap with the KN 
scattering states and the diquark-diquark field expected 
to have large enough overlap with the pentaquark state. 
To distinguish a single particle state from a scattering 
state the tools that we have at our disposal are the vol¬ 
ume dependence of the energy and spectral weights of the 
state provided we can isolate it with sufficient accuracy. 
The energy of a scattering state on a lattice, with the ex¬ 
ception of s-wave scattering states, is volume dependent. 
This is because, in the center of mass frame, the two par¬ 
ticles have non-zero relative momentum, which on a finite 
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lattice, depends on the spatial size of the lattice. This 
means that if a scattering state with non-zero momentum 
dominates the correlator then its energy will depend on 
the spatial length of the lattice. The other criterion that 
can be used is the scaling of the spectral weights with 
the spatial volume. For a two-particle scattering state 
the spectral weights are expected to be inversely propor¬ 
tional to the spatial volume whereas for a resonance there 
should be no volume dependence. The spectral decompo¬ 
sition of the correlator was given by Liischer who first 
noted that a weakly interacting two-particle state well 
below resonance energy should contribute to the correla¬ 
tor with an amplitude that scales inversely proportional 
to the volume. The scaling of the spectral weights as a 
probe for a two-particle resonance scattering state was 
first used in a lattice study of the pentaquark system by 
Mathur et al. Ea In order to check how reliably we 
can extract the scaling of spectral weights in practice we 
study the two pion system in the isospin 1=2 channel 
for which no low energy resonances are present. This 
study is done on three lattices of size 16^ x 32, 24^ x 32 
and 32^ x 64 at /3 = 6.0, which are the same lattices 
used in the study of the pentaquark system. We consider 
four different two pion interpolating fields and determine 
accurately the two lowest energy eigenstates. What we 
find is that the value of spectral weights extracted sta¬ 
bilizes only at large enough time separations and only 
then the ratio of spectral weights reaches the expected 
value. Therefore to check for volume dependence one 
requires very accurate data so that the extracted ratio 
of spectral weights is precise enough. For the two-pion 
system, given our statistics, this can only be achieved 
for the lowest eigenstate. How large contributions from 
higher relative momentum states are can be assessed by 
explicitly projecting to zero relative momentum. This 
is done for the lattice of size 16^ x 32 and the results 
obtained are compared to those without explicit projec¬ 
tion to zero relative momentum. We note that such a 
comparison has not been carried out in previous lattice 
studies of the pentaquark system. We confirm that di- 
agonalization of the correlation matrix without explicit 
projection yields a lowest energy eigenvalue that contains 
contributions from states with higher momenta and only 
becomes a pure s-wave scattering state at large time sep¬ 
arations. This means that in order to obtain the lowest 
scattering state one must allow for larger time separa¬ 
tion as the spatial extent of the lattice increases. As we 
will demonstrate for a lattice of spatial size of 3 fm one 
has to go beyond a time interval of 30a ^ 15GeH“^ to 
obtain the lowest scattering state. This is why we use 
Dirichlet boundary conditions (b. c.) in the temporal 
direction. In the study of the two-pion system we use 
the heaviest pion mass considered in this work, namely 
Ki = 0.153, to have the smallest statistical errors. We 
compare the scaling behavior of spectral weights in the 
two-pion and pentaquark system using the same value of 
Ki. Our main conclusion is that within the accuracy that 
the spectral weights can be determined from our data we 


can not exclude a pentaquark resonance. In the negative 
parity channel within our variational basis we only ob¬ 
tain the lowest eigenstate accurately enough to be able 
to perform the scaling analysis of the spectral weights. 
Because of this we are unable to draw a definite conclu¬ 
sion regarding the existence of the 0+. For the positive 
parity channel we obtain two eigenstates very close in en¬ 
ergy which at the two heaviest light quark masses have 
an energy gap of about 100 MeV independent of the spa¬ 
tial volume. As the quark mass gets lighter it becomes 
more difficult to obtain an accurate determination of the 
energy gap. Although the energy gap that we find is of 
the right order of magnitude the mass for this state is 
high: Measuring the mass in the positive parity channel 
for five values of light quark masses on our largest lattice 
using a smeared diquark-diquark interpolating field for 
the source and a local one for the sink and extrapolating 
to the chiral limit we obtain a value of 1.65 ± 0.09 times 
the mass of the non-interacting KN state. This value is 
too high to be identified with the ©''■(1540). The values 
of the light quark mass that we use in carrying out the 
extrapolation correspond to pion masses in the range of 
about 900-420 MeV. 


II. LATTICE TECHNIQUES 

In this work we consider two interpolating fields mo¬ 
tivated by recent proposals on the possible structure of 
the ©■'■ state. The first is based on the idea of diquark 
formation m and the other on a diquark-triquark struc¬ 
ture 0. Both have been used in previous lattice studies 
but a correlation matrix analysis in the manner consid¬ 
ered here has not been presented. We will refer to the 
first as diquark-diquark interpolating field. It is given by 

Jdd = Cs^ {<Cdf) (u^Cj^dh) (1) 

where C = 7072 is the charge conjugation operator. It 
was first proposed and used in the study of the pen¬ 
taquark system on the lattice by Sasaki [2. After per¬ 
forming the antisymmetric tensor contraction, we gen¬ 
eralize the resulting expression to obtain an isospin 1=0 
and 1=1 interpolating fields: 

Jdd = (ulC^^db) [{u^Cde) T {u^Cde)] CsJ . 

( 2 ) 

The minus sign corresponds to isospin 1=0 and the plus 
sign to isospin 1=1 respectively. The second field, which 
we will refer to as the KN-interpolating field is given by 

Jnk = (ulCjbdb) [ucisjbd) T ddsjbu)] , (3) 

where the minus sign corresponds to the isoscalar and 
the plus to the isovector. Although we have used both 
the isoscalar and the isovector interpolating fields, in this 
work we will only discuss the results in the isospin zero 
channel. 
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As presented, the interpolating fields have opposite 
parities with Jdd having positive parity and Jkn neg¬ 
ative. The parity of these fields can be flipped by mul¬ 
tiplying with 75 . The two point correlator can couple 
both to positive and negative parity states and the prop¬ 
agators for the positive and negative parity interpolating 
fields are given, respectively, by 

G+(f) = (4) 

G-(t) = 

g(t) = e{t)C+e-"^^*+ e{-t)C-e^~^ , 

where to"*" (m“) is the mass of the positive (negative) 
parity state. Employing Dirichlet boundary conditions 
means that only the terms with t > 0 contribute and 
therefore the positive and negative parity states corre¬ 
spond to the upper and lower Dirac components of 0+ (t) 
or the lower and upper components of G_ (t) respectively. 
For the two smaller volumes we use only Dirichlet b.c. 
and local interpolating fields since the main purpose of 
the evaluation on the smaller lattices is to examine the 
volume dependence of the spectral weights in the cor¬ 
relators. In order to perform our volume studies we 
also use Dirichlet b.c. for the local interpolating field 
on our largest volume. However the 32^ x 64 lattice has 
a time extent large enough so that backwards moving 
quarks are suppressed and anti-periodic boundary condi¬ 
tions can also be used without affecting the identification 
of the parity. Therefore on the large volume we opt for 
anti-periodic b.c. when using smeared sources. It is well 
known that smearing improves the overlap of the interpo¬ 
lating field with the ground state since it produces an ex¬ 
tended source of size typical to that of phy sical hadrons. 
We perform gauge invariant smearing 1241 by replacing 
a local quark field u(x) appearing in the interpolating 
fields by a smeared one, 'u(x), obtained by 

= ■ (5) 

y 

The gauge invariant smearing function $(x, y;t/(f)) is 
given by 

$(x, y; [/(f)) = (1 -h aff)^(x, y; U (t)) ( 6 ) 

where the hopping matrix H is defined by 

3 

H{x,y;U{t)) \Uj{x,t)6^ y_-. + u]{x - 

/=i 

(7) 

The parameters a = 4 and n = 50 are optimized to 
approximately reproduce the root mean square radius of 
the nucleon in the quenched theory at [3 = 6.0 where all 
the computations are carried out. 

Having M interpolating fields we can construct an M x 
M mass correlation matrix 

= J d^x< 0|Ji(x,f)j/(0,0)|0 > (8) 


where in the case of the pentaquark system the indices 
i,j = DD or KN. For non-vanishing results we must 
use interpolating fields of the same parity. This can be 
easily achieved by multiplying for example Jdd with 75 . 
Our variational analysis is performed in two ways: 

1. We solve the generalized eigenvalne equation 

C(fl)u„(fl) = A„(fi,fo)C(fo)'Cn(fl) • (9) 

For a large time separation ti — fp the eigenvalues are 
given by A„(fi,fo) = exp(-F;„(fi -fp)), n = 0 , 1 , 

1, yielding the energies En of the M lowest states with 
total momentum zero. The energy of the eigenstate 
is usually poorly determined since it has contributions 
from all the higher excited states, except when our vari¬ 
ational basis contains interpolating fields with a sizable 
overlap with all M lowest states. In this analysis we take 
fp/a = 3 — 6 with a the lattice spacing and check that 
the values that we find for En do not change as we vary 
to- 

2. We first diagonalize the correlation matrix C(fp) 

C{to)Vn{to) = A„(fo)w„(fp) (10) 

to determine the eigenvectors v{to) taking to/a = 1. We 
use these eigenvectors to project to the space spanned by 
the N largest eigenvalues A„(fp) 

Cj^it) = {vi,C{t)vj), i,j = 0,..,N-l (11) 

and solve the generalized eigenvector equation given in 
Eq.|21but for the projected correlation matrix C[j instead 
of C{t). We denote the resulting eigenvalues by A„. The 
eigenvectors, Vn, that we find determine the best linear 
combination Vn-Jn that has maximum overlap with 
the N lowest eigenstates. We will refer to this linear 
combination as the optimal interpolating field, JopUmai- 
The optimal correlator can be obtained by projecting 
C^: {yi,C^Vj) which yields the same energies as those 
extracted from the eigenvalues A„. We check that the 
eigenvalues, A„ that we obtain by diagonalizing are 
in agreement with A„. 

For the pentaquark system, in addition, to Cdd;KN 
constructed using the local interpolating fields J^dd and 
Jkn as a basis, we also consider Cdd-kn constructed 
using smeared fields Jdd and Jkn for the source but 
keeping the sink local and vice versa. The eigenvalues, 
A„, extracted from this correlation matrix should yield, 
for large enough time separations, the same energies as 
those extracted from A„. We check for consistency at the 
two heaviest light quark masses, namely ki = 0.153 and 
m = 0.155. 

The contractions needed for the computation of the 
pentaquark matrix elements are optimized by doing all 
the Dirac contractions explicitly before summing over 
color. For example, using the Jdd interpolating field 
as source and sink, the correlator is given by 
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CDD-,DD{t) = ^ 

X 


Tr (c-f 5 D‘^'<={x)Cj 5 U’’'^^{x)'^ Tr (cD^'f {x)CU<^'^(x)'^ 

- Tr (c-f 5 D‘^'^{x)CU’’''’^{x)'^ Tr (cd/'°(x)C' 75 C/“'“ ^(x)) 
+ Tr (c'75i:>^'"(x)C'75?7“'“'^(x)C'D^'^(x)C'i7^'^^(x)) 

- Tr (c'75i:>^'^(x)C?7“'“^(x)C'i:)^'^(x)C75C/'’"’^(x)) 



/ 

a ^ 




( 12 ) 


where S{x),U{x) and D{x) denote the strange, up and 
down quark propagators, Latin letters denote color in¬ 
dices, Greek letters Dirac indices and the transpose acts 
only on Dirac indices. In each of the terms we perform 
the Dirac algebra separately and store the result in an 
array e.g. in the first term having two traces we con¬ 
struct the traces and then we multiply them and sum 
over color indices. This reduces the time needed by more 
than an order of magnitude to about twice the time one 
needs for the nucleon correlator. Such grouping of terms 
is common in many lattice studies of matrix elements as 
for example in Refs. dllllll. 


III. RESULTS 

The lattices and values of the hopping parameter, ki, 
that we use for the light quarks are listed in Table ^ 
where we also give the mass of the pion, the kaon and the 
nucleon at each value of ki as well as the value obtained 
by linear extrapolation to the chiral limit using the form 

mjj = a + bml . (13) 


TABLE I: In the first column we give the lattice volume and 
in the second column the value of the hopping parameter, ki , 
for light quarks. In the third, fourth, fifth and sixth columns 
we give the pion mass, the kaon mass, the nucleon mass and 
the N* mass in lattice units. For ki = 0.1558 and 0.1562 
the mass of N* is not given since no clear plateaus could be 
identified. 


volume 

Ki 


rriK 

rriN 

mjv . 

local 

source and sink using 

Dirichlet b.c. and Jdd, Jkn 

16® X 32 

0.153 

0.422(2) 

0.363(2) 

0.787(8) 

1.039(56) 


0.155 

0.295(3) 

0.295(3) 

0.623(14) 

0.895(86) 

24® X 32 

0.153 

0.426(3) 

0.369(2) 

0.790(10) 

1.026(42) 


0.155 

0.302(3) 

0.302(2) 

0.644(12) 

0.877(57) 

32® X 64 

0.153 

0.420(2) 

0.361(2) 

0.788(6) 

0.987(33) 


0.155 

0.294(2) 

0.294(2) 

0.633(7) 

0.849(46) 

smeared 

source and local sink using anti-periodic b.c. and Job 

32® X 64 

0.153 

0.418(2) 

0.358(2) 

0.778(6) 

0.990(30) 


0.1550 

0.292(2) 

0.292(2) 

0.625(7) 

0.802(44) 


0.1554 

0.262(2) 

0.277(2) 

0.590(9) 

0.748(49) 


0.1558 

0.229(2) 

0.262(2) 

0.553(9) 



0.1562 

0.192(2) 

0.246(2) 

0.513(10) 



Kc =0.1571 

0. 

0.207(5) 

0.420(10) 


smeared source and sink using anti-periodic b.c. and Jdd 

32® X 64 

0.153 

0.418(2) 

0.360(3) 

0.789(9) 

0.973(39) 


In Eq. El rriH denotes the mass of any of the hadrons 
considered in this work. In Table H] we also include the 
mass of the negative parity baryon N* at the values of 
Ki where we could identify a plateau. For the strange 
quark we take Ks = 0.155, which produces a mass for the 
(p meson of 0.421(3) in lattice units. Using the mass of 
the nucleon in the chiral limit given in Table Q] we find 
that the ratio of the mass of the (p meson to the mass 
of the nucleon at the chiral limit is = 1.002 ± 

0.025, which is very close to the physical ratio of 1.087 
verifying that the value of Kg = 0.155 used is very close 
(within 10%) to the physical strange quark mass. We 
analyzed 202 configurations for the 16^ x 32 lattice 
and 100 configurations for each of the 24^ x 32 and 32^ x 
64 lattices. 


A. The nucleon sector 

In order to assess our methods of analysis we first ex¬ 
amine the results obtained for the nucleon using the stan¬ 
dard interpolating field 

Jn{x) = e“'’°(Mf(x) (775^6 (x))'Uc(x), (14) 

which can be local as given in Eq. El or constructed 
from smeared fields u and d instead of u and d. We 
denote the latter by In evaluating the nucleon cor¬ 
relator we can use Jn for the source and sink or use 
for the source and 77jv for the sink. In general, denot¬ 
ing with H the appropriate hadronic state, the correla¬ 
tor C{t) = f d^x < > computed using 
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local interpolating fields for the source and the sink is 
referred to as the local-local correlator, whereas the cor¬ 
relator C{t) computed using a local (smeared) field for 
the source and smeared (local) for the sink is referred to 
as local-smeared correlator. Finally the correlator, C'(t), 
computed using smeared fields for both the source and 
the sink is referred to as smeared-smeared correlator. In 
Fig.n we show results for the nucleon effective mass de¬ 
fined by 

meff (t) = - log (15) 

on the lattice of size 32^ x 64 with Dirichlet boundary 
conditions. On same figure we also show results obtained 
from C{t) with anti-periodic boundary conditions in the 
time direction. As can be seen smearing improves the 
overlap with the nucleon ground state resulting in an 
earlier plateau. All the errors shown on our figures and 
given in the tables are statistical and they are determined 
using a jackknife analysis. 

We fit the effective mass assuming two state dominance 
to the form 

1 -|- 

meff(t) = mo - log (16) 

where mo is the mass of the ground state, Am = mi — mo 
is the mass gap between the ground state and the first 
excited state and z is the ratio of the overlap of the in¬ 
terpolating field with the first excited state as compared 
to the ground state. For the nucleon sector mo gives the 
mass of the nucleon and Am the mass gap between the 
nucleon and the N(1440) state. The values of uiq and Am 
determined from fitting the effective mass derived either 
from C{t) or C(t) are in agreement. The errors on the 
mass gap Am however are large not allowing an accurate 
determination of the mass of the excited state. The fact 
that the ground state overlap is enhanced as compared 
to that of the first excited state when smearing is used, 
is confirmed by the decrease in the value of z by a factor 
of about twenty. Smearing both source and sink tends to 
suppress even further the excited state contributions and 
decorrelates the results at successive time slices. The 
mass extracted by fitting the smeared-smeared correla¬ 
tor, C{t), to a single exponential is mjv = 0.789(9) con¬ 
sistent with that extracted from the local-smeared corre¬ 
lator. However, as can be seen in Fig. ^ the statistical 
error on the smeared-smeared results at each time slice is 
larger as compared to the errors of the local-smeared re¬ 
sults due to the gauge noise introduced at the sink. Such 
errors will make it difficult to distinguish two closely ly¬ 
ing states and therefore in the case of the pentaquark 
system we will mostly use local-local and local-smeared 
correlators. 

An important test that distinguishes particle states 
from scattering states is the scaling of the spectral 
weights of local correlators with the spatial volume of the 
lattice [Hill. Expanding the correlator computed on a 



FIG. 1: The upper graph shows the effective mass for the nu¬ 
cleon (positive parity) and the lower graph for the N* (neg¬ 
ative parity) at k; = 0.153. The crosses show the results 
using the local-local correlator, C{t), and Dirichlet b.c. and 
the open triangles using the local-smeared correlator, C{t) 
and anti-periodic b.c. shifted to the right for clarity. In the 
case of the nucleon we also show results obtained using the 
smeared-smeared correlator C{t) (asterisks) shifted to the left 
for clarity. The dotted and dashed lines are fits to the effective 
masses obtained using C{t) and C{t) respectively assuming 
two state dominance. 

lattice of spatial size L in terms of states with the same 
quantum numbers as the interpolating field one obtains 

CL{t) = Y.W£e-^-^ . (17) 

a 

For a single particle state the weights are of order one 
whereas for a scattering state consisting of two weakly 
interacting particles well below the resonance are of or¬ 
der 1/L^. Furthermore, for a single particle state, 
determines the mass of the state and should be volume 
independent for large enough volumes. For a scattering 
state that is not an s-wave, on the other hand, is vol¬ 
ume dependent since it involves the relative momentum 
of the two particles, a fact which on the lattice means 
that each particle carries equal momentum in units of 
27r/L. In the case of the nucleon we know that we have 
a single particle state and we can check how reliably we 
can extract these weights. In Fig. [3 we plot the ratio of 
the local-local correlator CiQ{t) to C 24 (t) (computed on 
the lattices of size 16^ and 24^ respectively) for the nu¬ 
cleon and the N*. In a time interval where a single state 
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FIG. 2: We show the ratio /C 2 i{t) (crosses) and the 

ratio i?i 6;24 (open triangles with a small shift to the right for 
clarity) as a function of the time separation in lattice units. 
The upper graph shows results for the nucleon and the lower 
for N* at Ki = 0.153. 


dominates (plateau region), assuming its mass is volume 
independent, the ratio of correlators is equal to the ratio 
of weights Wi^jWiA- This ratio should be one for a single 
particle state. As can be seen the ratio of correlators is 
approximately one. In the same figure we also show the 
ratio 


R 


H 

Li :L2 




(18) 


where with H we denote the nucleon or the N* and we 
have taken Li = 16 and L 2 = 24. This ratio corrects 
small finite size effects on the mass and in the plateau 
region is again equal to ■ As can be seen 

indeed improves the signal, in particular for the nucleon, 
giving better agreement with unity, which is the expected 
result. 


B. The two pion sector 

Having tested our techniques for a single particle state 
we now verify that they can be applied for two particle 
scattering states. For this test we choose the two pion 
system in the isospin two channel where no low lying 
resonances are expected. We consider five interpolating 


fields: 


Ji^x) 


Jt{^) 


Jl{,x)J^{x), (x) = d{x)-f 5 u{x) 

J2{x)J2{x), J 2 = dix)-f5'you{x) 
J3ix)J^{x), = d{x)-f5e^-/f,u{x) 

3 

Jo{x)J^{x), = d{x)'yo'^'yiu{x) 

i=l 

Jf{x)Jf{x), Jf{x) = d{x)-f^u{x) (19) 


The first three are products of pion interpolating fields 
whereas the last two are products of rho-type interpo¬ 
lating fields. As we already pointed out the energy of a 
scattering state of two hadrons hi and /12 depends on the 
spatial size L of the lattice and is given by 

Eh^h2 = j (t) ’ = 0’ ■ 

( 20 ) 

where we have suppressed the L index on the energy. In 
Fig.El we show the effective mass for a single pion using 
the interpolating fields Ji{x), J^ix)^ and J^[x) for our 
three lattices. As expected Ji{x), routinely used in lat¬ 
tice calculations, and its variant J 2 {x) have the largest 
overlap with the pion. All three yield consistent results 
for the pion mass. For the two pion system the inter¬ 
polating field Ji^{x) is very noisy as compared to the 
other four and we do not include it in the figure. We 
obtain the best overlap with the two pion ground state 
when using Ji^{x) with Jl'^(x) and J^'^{x) being the 
next best. Again all interpolating fields produce consis¬ 
tent two pion energies for large enough time separation. 
For lattices 16^ x 32 and 32^ x 64 the time extent is suf¬ 
ficient to have a large plateau region after suppression 
of excited state contributions. For the 24^ x 32 lattice 
this occurs for time separations t/a > 22 limiting the fit 
range to a few time slices. 

We perform a variational analysis using the interpolat¬ 
ing fields defined in Eq. excluding , which means 
that in Eq. El we take M = 4. To check we also compute 
the eigenvalues by the second method of analysis namely 
by projecting using Eq. CH to an A^ X matrix where 
we take N = 3. What we find is that we obtain two 
distinct eigenstates for which we also construct the best 
interpolating field using the eigenvectors V„. For this sys¬ 
tem had we used only interpolating fields l7i and J 4 the 
resulting eigenvalues would be consistent with the ones 
obtained within the larger variational set. This is an im¬ 
portant observation since it demonstrates that using two 
interpolating fields which have very good overlap with 
two low lying states we can accurately determine them 
by diagonalizing the 2x2 correlation matrix. The effec¬ 
tive mass determined for the two eigenvalues A^, is shown 
in Fig. 21 for our three lattices. In all cases the effective 
mass extracted from the lowest energy eigenvalue (largest 
A) is in agreement with that extracted from the correla¬ 
tor that uses the interpolating field . As the spatial 
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FIG. 3: The effective mass at ki = 0.153 as a function of t/a 
for the lattice of size 16® x 32 (top), 24® x 32 (middle) and 
32® X 64 (bottom). The set results with the lower value are 
obtained using one pion interpolating fields, Jx (crosses), J 2 
(asterisks) and (open triangles). The set with the higher 
value correspond to the two pion interpolating fields, 
(crosses), (asterisks), (open triangles), and (cir¬ 
cles). The dotted lines are the two lowest two-pion scattering 
states i? 27 r and 


FIG. 4: The effective mass at ki = 0.153 extracted from the 
correlator using the interpolating field (crosses) and from 
the two largest eigenvalues Aq (open triangles) and Ai (aster¬ 
isks) for the lattice of size 16® x 32 (upper graph), 24® x 32 
(middle graph) and 32® x 64 (lower graph). The dotted lines 
are the two lowest energy two-pion scattering states and 
E 2 Tr- The dashed line is the lowest energy two-rho scattering 
state E 2 p. 


extension of the lattice increases a larger time separa¬ 
tion is needed for the effective mass to reach the energy 
of the s-wave scattering state Since in this system 

the low lying states are two particle scattering states this 
clearly demonstrates that the two lowest eigenvalues do 
not correspond to a given value of the relative momen¬ 
tum carried by each particle but they are an admixture 
of states having non zero relative momentum that mix 
with the state with zero relative momentum. The energy 
difference, if" — between different relative momen¬ 

tum states decreases like 1/L^ and therefore a larger time 
interval is needed in order to reach the ground state as 
L increases. This can be seen in Fig. 0] where for the 
16® lattice the plateau for the ground state is reached 
for t/a > 13 whereas for the 24® we need t/a > 23 and 
for the 32® t/a > 30. These values determine the small¬ 
est value of the lower time range to be used in fitting 
the effective masses to a constant. If one would like to 
reduce the lower fit range then fitting to a two exponen¬ 
tial form using Eg. 1161 is essential. A similar behavior is 


also observed for the second eigenstate: On the 16® the 
two-rho scattering state E^p is the second lowest energy 
state since E^^ is higher and a clear plateau is seen for 
t/a > 17. Going to the larger 24® lattice E^^^ decreases 
becoming about equal to E^p with a plateau that sets 
in for t/a ^ 25 resulting in a fit range limited to three 
points. This explains why the energy extracted is higher 
on this lattice than 2mp as can be seen in Fig. [S] where 
we plot the mass extracted from fitting the effective mass 
in the available plateau region versus the spatial lattice 
size. Finally, despite the fact that on the 32® lattice 
E 2 t^ < El/p the two-rho state dominates the time depen¬ 
dence of the second eigenvalue for 30 < t/a < 40. For 
t/a > 40 the effective mass becomes consistent with E/.^ 
but in this time range the statistical errors have become 
much larger. Therefore fitting the effective mass in a 
time range 30 — 40 or extending the upper fit range af¬ 
fects very little the resulting value thereby obtaining the 
rather accurate measurement of E/p displayed in Fig.0 
Note that corrections due to interactions between the 
























FIG. 5: The mass at k; = 0.153 extracted from fitting 
the plateau of the effective mass of the lowest energy eigen¬ 
value (filled triangles) and the second higher energy eigen¬ 
value (Hlled circles shifted to right for clarity) as a function 
of the spatial length of the lattice L. Also shown are results 
for twice the mass extracted from the pion correlator (crosses) 
and rho correlator (asterisks). The dotted line is the two-pion 
scattering energy E 2 „. 


two pions cannot be seen on the scale of this figure as¬ 
suming that we can use the physical scattering length in 
Luscher’s result [2^. It must be noted that on the large 
lattice fitting in the plateau-like region in the time range 
15 < t/a < 30 would yield an incorrect higher value for 
the energy. Therefore this demonstrates that for lattices 
of spatial extent of about 3 fm we need a time separation 
of at least 40a to obtain correctly the energy of the two- 
rho state making the use of Dirichlet b.c. essential. To 
verify that the reason the plateau region starts at relative 
large time separations is due, to a large extent, to con¬ 
tamination of states with higher relative momentum, we 
perform on the small lattice an explicit momentum pro¬ 
jection to two particles each carrying zero momentum. 
This is done by evaluating the correlation matrix 

= E < 0|j/(x,t)j/(y,t)4^t(0)4^^(0)|0 > , 

x.y 

( 21 ) 

where s, s' denotes the pion or the rho. 

In Fig.Elwe compare the eigenvalues obtained when we 
carry out the zero momentum projection for each particle 
to our previous (unprojected) results. As can be seen 
for both eigenstates the plateau region starts at earlier 
time separations. In fact for t/a > 6 the projected two 
particle correlator is the same as the product of the single 
particle correlators. For large enough time separations 
the unprojected eigenvalues approach the correct values 
E 2 ^ and i? 2 p. The mass gap E^p — E^^ ~ 320 MeV 
at Ki = 0.153. The mass gap between the 0''‘ and the 
KN s-wave scattering state is not known at this value of 
Ki- As will be discussed below, the mass gap between 
our candidate pentaquark resonance and increases 



FIG. 6: The upper graph shows the effective mass for the low¬ 
est energy eigenstates and the lower graph the effective mass 
for the second higher energy eigenstate. In both graph we 
show data with projection to zero relative momentum using 
Eq. 1211 (asterisks shifted to the right for clarity) and with¬ 
out explicit projection (open triangles) for the lattice of size 
16® X 32 at Ki = 0.153. We also show twice the effective mass 
extracted from a single pion (crosses in upper graph ) and 
rho correlator (circles in lower graph). The dotted line is the 
two-pion scattering state E^-w- 


with the quark mass to about 170 MeV at ki = 0.153 in 
the negative parity channel. This estimated smaller gap 
can make the study of the pentaquark system harder. 

Having identified the two lowest eigenstates within our 
variational basis we study the spectral weights of these 
states. In Fig. [ 7 | we show the ratio i?i 6:24 for the pion 
when using the interpolating field (x) and for two pi¬ 
ons when using as well as ^optimal- This ratio 

is unity for the pion as expected whereas for the two- 
pion system it increases approaching the expected ratio 
of 3.4 for ti/a > 25 when the s-wave two pion scattering 
state dominates. For the lattice of spatial extent 32 this 
happens for t/a > 30 and therefore i?i 6:32 stays close to 
unity up to about t/a = 30 which is the largest time 
separation for which it can be constructed. Instead of 
7 ?i 6:32 we can extract the spectral weights by fitting the 
correlator to a single exponential or to a sum of two ex¬ 
ponentials. This allows to take into account information 
from the full time extent of the lattice. We fix the upper 







9 



t/a 


FIG. 7: The ratio i?i 6:24 at ki — 0.153 for one and two pion 
states using interpolating field J\(x) (crosses) and {x) 
(asterisks) as a function of tja. The ratio i?i 6;24 extracted 
from correlators using the optimal combination Joptimai{x) 
for the lowest eigenstate of the two pion system is shown by 
the open triangles shifted to the right for clarity. The dotted 
lines show the expected value of the ratio for a single particle 
state and two-particle scattering state. 


fit range to 26 in lattice units for the lattices of temporal 
extent Nt = 32 and to 56 for the lattice with Nt = 64. 
We show in Fig.|Slthe ratio of spectral weights 
and VF 16 /W 32 for the lowest state. What can be seen 
is that both ratios extracted using a single exponential 
fit increase approaching the expected value for a scat¬ 
tering state. Fitting to a sum of two exponentials we 
obtain ratios that approach the expected value at much 
smaller time separations albeit with larger errors. How¬ 
ever if instead of 56 we take an upper fit range of 26 for 
the large lattice either using single or double exponential 
fits the ratio FF 16 /W 32 stays close to unity leading to the 
wrong conclusion. We show the corresponding ratios for 
the second eigenstate in Fig. 0 using a single exponen¬ 
tial. Even though the errors on the correlators are small 
the error on the ratio of spectral weights is too large to 
lead to a definite conclusion. Had we used a sum of two 
exponentials the errors would be even larger. 


C. The pentaquark sector 

The first experimental indication for the existence of a 
pentaquark state came when the LEPS collaboration P| 
detected a resonance with mass 1540 ±10 MeV and width 
less than 25 MeV in accord with the predictions of the 
chiral soliton model This small width is surprising 
since, lying about 100 MeV above the KN threshold, it is 
expected to readily decay to KN. Since the observation of 
the first signal for 0 ''" was reported, several other experi¬ 
ments were carried out, some confirming its existence 0 
and others not ^ . 

Whether lattice QCD supports the existence of a pen- 
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FIG. 8: The ratio of spectral weights 'u;i6/ui24 (upper graph) 
and wie/W 32 (lower graph) as a function of the lower fit range 
ti/a at Ki — 0.153. Asterisks denote results extracted from 
the correlators with {x), open triangles and filled circles 
from correlators with JopUmai (x) fitted to a single exponential 
to a sum of two exponentials respectively. For comparison we 
also show these ratios for the single pion state using J\(x) 
(crosses). The upper fit range is fixed to 26 for the 16® x 32 
lattice and to 56 for the 32® x 64 lattice. On the lower graph 
we show with the filled triangles the ratio of spectral weights 
W 32 when the upper fit range for the large lattice is fixed 
to 26. The dotted lines show the expected value of the ratios 
for a single particle state and a two-particle scattering state. 


taquark state is also unclear [^. The study of reso¬ 
nances in lattice QCD is in its infancy and therefore it 
comes with no surprise that identifying the on the 
lattice turned out to be a difficult task. In the two pion 
system scaling of the spectral weights W with the spatial 
volume was observed for large time separations. We will 
use the largest available time separations in applying the 
same analysis in the pentaquark system. Clearly whether 
this criterion can be used in practice will depend on the 
size of the statistical errors. 

In Figs.irni andllllwe show the effective mass for both 
the negative and positive parity states evaluated on the 
three lattices at ki = 0.153 for isospin zero using the lo¬ 
cal interpolating fields Jdd and Jkn- What can be seen 
from these figures is that both interpolating fields yield 
consistent plateaus in both channels with the KN inter- 
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FIG. 9: The ratio of spectral weights w\Q/w 2 i extracted from 
the correlators with interpolating held JopUmaiix) (asterisks) 
for the hrst excited state as a function of the lower ht range 
ii/a. The corresponding ratio UI 16 /UI 32 is shown with the 
open triangles. The dotted lines show the expected value of 
the ratio for a single particle state and two-particle scattering 
state. The crosses show for comparison the ratio of spectral 
weights for one pion state extracted from the correlator using 
interpolating held J^{x). 


polating field yielding smaller errors in the negative par¬ 
ity channel than the diquark-diquark interpolating field 
whereas in the positive parity channel the opposite is 
true. We also see that the errors even at this heavy pion 
mass are large especially for the positive parity channel 
where the large time behavior is not well determined. 
This means that extracting the spectral weights for the 
positive parity channel accurately enough for large time 
separations to test scaling will not be possible. There¬ 
fore for the positive parity channel we can only look for 
volume dependence in the extracted energy of the states. 
This energy can be determined either by fitting the effec¬ 
tive mass in the plateau range to a constant or, allowing 
for two states, to the form given in Eq. El The latter is 
especially important when using local-local correlators, 
where excited states have a large contribution up to time 
separation of about 14a. Since the aim is to detect a res¬ 
onance, which is higher than the scattering KN state, it 
is crucial to be able to fit in a time range where it is still 
the dominant state. This is expected to be so for time 
intervals less than about t/a ~ 20 or t ~ 10 GeV“^. The 
values obtained for the mass, fitting within this range us¬ 
ing either Jdd or Jkn, are given in Table ITU for all three 
lattices. From these values, we can conclude that either 
using the diquark-diquark or the KN local interpolating 
field the masses that we find are within each other’s er¬ 
ror. The systematic volume effect of a decrease in the 
mass with the increase of the lattice size, characteristic 
of a scattering state, is not observed in the positive parity 
channel. 

The standard approach to suppress excited state con¬ 
tributions improving the effective mass plateaus is to use 



FIG. 10: The effective mass for the pentaquark in the negative 
parity channel for lattices of size 16^ x 32 (upper graph), 24® x 
32 (middle graph) and 32® x 64 (lower graph) at ki = 0.153. 
The crosses show the results obtained from correlators using 
J^DD and the open triangles using Jkn- The dotted lines 
show mjv -I- niK and 


smearing. We use Wuppertal smeared interpolating fields 
on the large lattice imposing antiperiodic boundary con¬ 
ditions in the temporal direction. As can be seen in 
Fig. 1121 where we show results for ki = 0.153, smearing 
the source but keeping the sink local suppresses excited 
state contributions. Note that, within our statistical er¬ 
rors, the effective mass with smeared source on half the 
lattice is in agreement with the results from the local cor¬ 
relator at larger times leading to the same plateau values. 
This justifies the use of antiperiodic b.c. As we already 
pointed out if the sink is also smeared the gauge noise 
increases making identification of two close-by states im¬ 
possible. Therefore for the rest of the results we will 
apply smearing only at the source. 

Having confirmed that both local and smeared inter¬ 
polating fields as well as the KN and diquark-diquark 
interpolating fields lead to consistent values for the mass 
we perform the variational analysis described in section 
II. As in the two-pion system the eigenvalues will still 
be contaminated by scattering states with non-zero rela¬ 
tive momentum. This is particularly severe on the large 
lattice where the scattering KN states are very close in 
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FIG. 11: The effective mass for the pentaquark in the positive 
parity channel at ki = 0.153. The dotted lines show s-nd 
i?x]v fl'iid the dashed line mx + mjv*. The rest of the notation 
is the same as that of Fig. ITUl 


energy as can be seen in Figs. ^1 and m If within our 
variational basis we have interpolating fields that have 
good overlap with the pentaquark resonance then the two 
energy eigenvalues should yield the KN scattering states 
and the 0'*". The choice of the diquark-diquark inter¬ 
polating field is based on the expectation that it should 
approximate well the structure of 0“''. However if this 
is not the case then we will not be able to obtain the 
resonance state as the second eigenvalue of the 2x2 cor¬ 
relation matrix. The lowest eigenstate will be accurately 
determined and, allowing for large enough time separa¬ 
tions, it should produce the s-wave KN scattering state in 
the negative channel and the p-wave KN scattering state 
in the positive channel. If there is an admixture of a reso¬ 
nance state then we expect the spectral weights to show a 
different behavior from that observed in the two-pion sys¬ 
tem where the only low lying states are scattering states. 
In the positive parity channel only KN scattering states 
with non-zero relative momentum contribute. Therefore 
the lowest energy is expected to decrease as the spatial 
volume of the lattice increases according to Eq. HOI Since 
we expect the errors on the spectral weights to be too 
large in this channel to allow us to test scaling the only 
option is to study the volume dependence of the energy 



negative parity channel and the lower graph for the positive 
parity channel at ki = 0.153 on the lattice of size 32® x64. The 
crosses show results obtained from the local-local correlator 
C{t), the open triangles results from the local-smeared corre¬ 
lator C{t) shifted to the right and the asterisks results from 
the smeared-smeared correlator C{t). The diquark-diquark 
interpolating field is used in all correlators. The dashed line 
is a two-exponential fit to the local-smeared results. The dot¬ 
ted line shows FIxat for the negative and for the positive 
parity channel. 


on the spatial length of the lattice. However in the pos¬ 
itive parity channel the s-wave KN* scattering state is 
allowed. The energy of this state is shown in Fig. CH On 
the small lattice the value of the smallest relative momen¬ 
tum 27r/L is such that, at ki = 0.153 using the masses 
given in Table Ql comes out very close to Ej^j^ 

whereas for the two larger lattices is higher. Therefore if 
the two lowest eigenstates are the KN* and KN scatter¬ 
ing states then we expect the energy gap — E^j^, \ 

while being almost zero on the small lattice to increase for 
the other two. Clearly the existence of the KN* scatter¬ 
ing state complicates the identification of a pentaquark 
resonance also in the positive parity channel and would 
require at least a 3 x 3 or even a 4 x 4 correlation matrix 
to allow us to resolve three different states. A very accu¬ 
rate determination of the spectral weights for the three 
lowest states will be needed in order to distinguish the 
KN and KN* scattering states from a resonance. This 
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FIG. 13: The upper graphs show the effective mass for the 
negative parity and the lower graphs for the positive parity 
channels at k; = 0.155 on the lattice of size 32® x 64. The 
squares on the left two graphs denote the results extracted 
from the correlators with the local interpolating field Jdd- 
The crosses on all four graphs show the effective mass, de¬ 
rived from the eigenvalue Ao(t) of Cdd-,kn- On the graphs 
on the right the open triangles denote results derived from 
Ai(f), whereas the circles and asterisks are results extracted 
from the eigenvalue Ao(f) and Ai(t) from the diagonalization 
of CddiKN respectively. The dotted lines show i?xAr for the 
negative parity channel and for the positive parity chan¬ 
nel. 


is beyond the scope of this work. Our working hypoth¬ 
esis is that if the 0+ is present then it should dominate 
the correlator in the appropriate time range and to have 
largest coupling to the KN scattering state. 

In Fig. 1131 we show effective masses determined from 
the two eigenvalues Ao(t) and Ai(t) of the correlation 
matrix C(t)uD-,KN for ki = 0.155 corresponding to a 
lighter quark mass than ki = 0.153 discussed up to now 
although the behavior is similar. We observe that, like 
in the case of the two-pion system, the effective mass 
derived from the lowest energy eigenvalue is consistent 
with that obtained from either interpolating field J^dd 
or Jkn- On the same figure we also show the effective 
mass obtained from the eigenvalues Ao(t) and Ai(t) of 
the correlation matrix C(t)]jo-KN, which involves local- 
smeared correlators of the two interpolating fields Jdd 
and Jkn and antiperiodic boundary conditions in the 
temporal direction. Again both local-local and local- 
smeared results yield the same plateaus with smearing 



FIG. 14: The energy gap between the two lowest energy eigen¬ 
states as a function of t/a for the pentaquark in the positive 
parity channel at m = 0.153 (upper graph) and ni — 0.155 
(lower graph) for lattices of size 16® x 32 (asterisks), 24® x 32 
(crosses) and 32® x 64 (open triangles). The dash line shows 
the plateau value determined on the lattice of size 16® x 32. 


suppressing excited state contributions. The main obser¬ 
vation is that the first excited state behaves differently 
in the positive and negative parity channels: In the neg¬ 
ative parity channel is very high in energy and poorly 
determined. This is puzzling since if the diquark-diquark 
operator has good overlap with the 0“'' then we would 
have expected it to be the first excited state of this cor¬ 
relation matrix with the ground state being the KN scat¬ 
tering state. In fact if we compare the overlap of Jdd 
with the ground state to that obtained with Jkn we find 
that it is about 50 times smaller. This raises questions as 
to how well a local diquark description approximates the 
structure of a pentaquark resonance. In the positive par¬ 
ity channel we obtain two distinct eigenstates very close 
together. From our study of the two pion system we know 
that the eigenvalues, in general, will not have a definite 
relative momentum but will be contaminated with states 
with higher relative momenta. Therefore we do not ex¬ 
pect that for the time separations that we can determine 
the gap to be able to relate it to — E]^^. What 
we can do is to check if this gap depends on the spatial 
size of the lattice. We evaluate this energy gap on our 
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FIG. 15: The upper graph shows the effective mass for the 
lowest energy eigenstate of the two-pion system as a function 
of t/a at Ki = 0.153. The middle graph shows the correspond¬ 
ing effective mass for the pentaquark system in the negative 
parity channel and the lower graph for the pentaquark in the 
positive parity channel for the lattices of size 16^ x 32 (aster¬ 
isks), 24® X 32 (crosses) and 32® x 64 (open triangles). 


three lattices at two values of ki = 0.153 and ki = 0.155 
and plot the results in Fig.^l Within our statistical ac¬ 
curacy the energy gap does not show any strong volume 
dependence having about the same values of 100 MeV at 
the two K values. Although this energy gap agrees with 
the energy difference between the and the KN scat¬ 
tering threshold the mass that we estimate below for the 
positive channel at the chiral limit is much larger than 
1540 MeV to be identified as the 0+. 

Fitting the effective mass derived from Ao(t) to Ea.ITCl 
we find at ki = 0.153 Eq = 1.228(14) for the negative 
parity and Eq = 1.438(60) for the positive parity chan¬ 
nels whereas at ki = 0.155 we find Eq = 1.014(20) and 
Eq = 1.259(116). These values agree both with those 
extracted from the local-smeared correlator C{t) using 
Jdd given in Table cni and from Ao{t) as demonstrated 
in Fig. Ha The fact that these different ways of con¬ 
structing the correlation matrices lead to the same value 
for the energy makes us confident that our results for the 
lowest state are robust. For the second state the only 
estimate that we can provide is that, at our two heaviest 



FIG. 16: The ratio i?i 6:24 versus t/a for the negative par¬ 
ity channel for the lowest energy eigenstate at k; = 0.153. 
Asterisks show results obtained with correlators using the 
diquark-diquark interpolating field Jdd and open triangles 
with ^optimal- For comparison we also show the ratio for the 
nucleon (crosses). The dotted lines show the expected value 
of the ratios for a single particle state and a two-particle scat¬ 
tering state. 


pion masses, is higher in energy by about 100 MeV in the 
positive parity channel. 

In Fig. E] we compare the effective masses for the low¬ 
est state of the pentaquark system and of the two-pion 
system. The effective masses in the two-pion system start 
off with higher values as the spatial volume increases be¬ 
cause scattering states with higher relative momentum 
have smaller energy gap at larger spatial volumes, re¬ 
quiring longer time to damp out. Only at large times 
the ground state dominates the correlator yielding the 
same plateau value on all lattice sizes. In the pentaquark 
system, on the other hand, the effective masses have the 
same value in the time interval 20 < t/a < 30. This is the 
time interval where a resonance is expected to dominate 
the correlator. However one has to bear in mind that the 
statistical errors are bigger than in the two-pion system 
so that any volume dependence maybe harder to detect. 
This is particularly true in the positive parity channel. 
Therefore we check if scaling of the spectral weights can 
yield information on the nature of these states. For the 
same reasons as explained for the two-pion system we 
only show the ratio of correlators i?i 6 : 24 - As can be seen 
in Fig. [TBl for the negative parity channel i?i 6:24 stays 
close to unity up to time separations t/a ~ 20. For 
t/a > 20 where, from our study of the two-pion sys¬ 
tem, we expect this ratio to start to deviate from unity 
the errors become large making it difficult to distinguish 
a single particle state from a scattering state. In the 
positive channel in this large time window the errors are 
even larger making this test not applicable. Instead of 
the ratio i?i 6:24 we plot the spectral weights determined 
on each volume making a direct comparison to the two- 
pion system. Again we do this only for the negative par- 
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FIG. 17: We show the spectral weights for the lowest eigen¬ 
state at Ki = 0.153 versus the lower htting range ti/a keeping 
the upper range fixed to 26 for the two smaller lattices and 
to 56 for the large lattice. The upper graph shows results for 
the two-pion system and the lower graph for the pentaquark 
in the negative parity channel for the lattices of size 16^ x 32 
(asterisks for a single exponential fit), 24^ x 32 (crosses for 
a single exponential fit) and 32® x 64 (open triangles for a 
single exponential fit). In the case of the two pion system 
we also show results for the weights determined from fitting 
the correlator to a sum of two exponentials for lattices of size 
16® X 32 (squares), 24® x 32 (circles) and 32® x 64 (filled tri¬ 
angles). The solid lines in the upper graph show the plateau 
values. 

ity channel where we have more accurate results. This 
comparison is shown in Fig. El as a function of the lower 
time range U/a used in the single and double exponential 
fits to the optimal correlators keeping the upper fit range 
fixed to 26 for the two smaller lattices and to 56 for the 
large lattice. For the two-pion system the weights are 
clearly volume dependent. On the lattice of size 16^ x 32 
the values we obtain, performing a single exponential fit, 
are independent of the lower fit range for ti/a > 10 and 
we can reliably extract the spectral weight of the lowest 
state. For the 24® x 32 lattice the convergence is seen 
later because suppression of higher momentum states re¬ 
quires larger time separations. For the 32® x 64 lattice 
the spectral weights are less well determined and can vary 
by a factor of about two depending which range of ti/a 


we use. However, within this uncertainty, the spectral 
weight ratio is closer to the values expected for a scat¬ 
tering state than to unity. Fitting to a sum of two expo¬ 
nentials the spectral weights become independent of the 
lower fit range for ti/a > 10 on all three lattices yielding 
results consistent to those obtained when fitting to a sin¬ 
gle exponential. This checks that indeed the value we find 
for the spectral weights does not depend on our fitting 
scheme. The volume dependence of spectral weights for 
the pentaquark system is different. The spectral weights 
extracted from a fit to a single exponential are vol¬ 
ume independent for 10 < ti/a < 26, which is the time 
range where we can make a comparison. This means that 
in this range we do not have a single scattering state. The 
energy gap between the two lowest KN scattering states 
for Ki = 0.153 is a{Ej^j^ — = 0.26 and 0.13 on 

the lattices of size 16® x 32 and 24® x 32 respectively as 
compared to the energy gap of a{2mp — 2m,n.) = 0.16 in 
the two-pion system. This means that, as compared to 
the two-pion system, the lowest KN scattering state is 
filtered out at smaller time separations on the smallest 
lattice and at ~ 20% larger times on the 24"^ x 32 lattice. 
Therefore one would have expected to see a clear volume 
dependence of the spectral weights on the two smaller 
lattices which is not observed. On the largest lattice the 
fact that the value of spectral weight decreases as ti/a in¬ 
creases beyond 30 is consistent with the expectation that 
at large time separations a KN scattering state should 
dominate. For the positive parity channel extending the 
upper fit range to large time separations yields statistical 
errors that are too large to make this test conclusive even 
at this very heavy pion mass. As the light quark mass 
decreases the errors on the determination of the spec¬ 
tral weights increase and therefore one would need much 
larger statistics to reliably extract the weights. Based on 
the volume study of the spectral weights at m = 0.153 
we can not exclude a pentaquark resonance. 

Within our variational basis the lowest energy eigen¬ 
value is the only one that we can determine accurately in 
the negative parity channel. It is shown to be in agree¬ 
ment with the mass extracted from the correlators using 
either local interpolating fields Jdd, Jkn and JopUmai 
or the smeared versions of these. Therefore to get an es¬ 
timate on the light quark mass dependence of the lowest 
state we use the local-smeared correlators C(t). Similarly 
for the positive parity channel the energy gap between the 
two lowest energy eigenvalues is about 100 MeV at both 
K = 0.153 and k = 0.155. Resolving these two states 
for smaller quark masses is even harder and therefore we 
opt to evaluate just the lowest energy which can again 
be obtained from just using Jdd- This will give us a 
rough idea of what the chiral limit of this state is. Since 
the local-smeared correlators are evaluated on our largest 
volume with periodic b.c. in the time direction the mass 
is extracted using the form given in Eg. 1 161 for which half 
the time separation is sufficient as demonstrated by the 
agreement of our results at ki = 0.153 and 0.155. 
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TABLE II: The first column gives the interpolating field, the second column gives the mass of the pentaquark extracted by 
fitting the effective mass to a constant in the plateau region and the third column gives the mass extracted by fitting to the 
form given in Eq. for the lattice of size 16® x 32. The corresponding quantities for the other two lattices are given in the 
next four columns. We give results at values of k; = 0.153 and 0.155 for both parity states. The /d.o.f for these fits is less 
than one. 



16® X 32 

24® X 32 

32® X 64 


K.I = 0.153 


Negative Parity 

operator 

am (1 exp) 

am (2 exp) 

am (1 exp) 

am (2 exp) 

am (1 exp) 

am (2 exp) 

Jdd 

1.172(14) 

1.174(14) 

1.236(14) 

1.240(15) 

1.237(12) 

1.234(13) 

•Jkn 

1.173(9) 

1.180(9) 

1.232(15) 

1.239(17) 

1.239(9) 

1.236(9) 

^optimal 

1.178(14) 

1.176(16) 

1.232(14) 

1.253(23) 

1.235(9) 

1.230(9) 


Positive Parity 

Jdd 

1.432(50) 

1.376(59) 

1.485(30) 

1.449(48) 

1.382(30) 

1.423(25) 

Jkn 

1.488(51) 

1.411(87) 

1.479(44) 

1.434(64) 

1.400(36) 

1.414(35) 

^optimal 

1.502(135) 

- 

1.448(52) 

1.465(167) 

1.462(46) 

1.444(68) 


K.I = 0.155 


Negative Parity 

operator 

am (1 exp) 

am (2 exp) 

am (1 exp) 

am (2 exp) 

am (1 exp) 

am (2 exp) 

Jdd 

0.977(44) 

0.938(23) 

1.014(19) 

1.013(23) 

1.031(15) 

1.034(15) 

Jkn 

0.947(15) 

0.937(17) 

0.998(20) 

1.002(26) 

1.023(10) 

1.020(11) 

i^optimal 

0.952(19) 

0.929(27) 

1.006(15) 

1.012(32) 

1.022(13) 

1.020(14) 


Positive Parity 

Jdd 

1.300(67) 

1.184(92) 

1.323(33) 

1.316(44) 

1.308(26) 

1.250(49) 

Jkn 

1.374(60) 

1.283(129) 

1.274(48) 

1.255(67) 

1.230(42) 

1.218(47) 

KToptimal 

1.457(89) 

- 

1.256(57) 

1.283(179) 

1.290(75) 

1.259(116) 


The results that we obtain for the mass from the local- 
smeared correlators for all the light quark masses that we 
have considered in this work are given in Table ITm We 
note that at ki = 0.155 diagonalization of the correlation 
matrix CoD-KNit) yields in the negative parity channel 
Eq = 1.002(15) for the ground state and = 1.6(2) 
for the first excited state. This value of Eq is in agree¬ 
ment with the value given in Table cni and the value of 
E^ is in agreement with 1.7(1) extracted from the local- 
smeared correlator via Eq. M The pentaquark masses 
at the chiral limit given in Table cni are obtained by lin¬ 
early extrapolating the results determined on the set of 
five Ki values using Eq. 1131 The corresponding results for 
the kaon and nucleon masses at the chiral limit are given 
in Tabled We show all chiral extrapolations in Fig. II81 
Note that the slope in the case of pentaquarks is steeper 
than for the nucleon or the kaon which means that the 
mass gap between pentaquarks and KN increases with 
the quark mass, giving at ki = 0.153 a gap of about 
170 MeV as mentioned earlier. From the values obtained 
at the chiral limit we can evaluate the ratios of the mass 
of the candidate pentaquark in the positive and nega¬ 
tive parity channels to the mass of the kaon-nucleon sys¬ 
tem, rriKN- The values that we find for these ratios are 
1.12 ± 0.04 for the negative parity and 1.65 ± 0.09 for 
the positive. Substituting the physical kaon and nucleon 


TABLE III: We give, in lattice units, the mass of the pen¬ 
taquark in the negative and positive parity channels, m~ and 
m^, determined using Eq. cni to fit the effective mass ex¬ 
tracted from the local-smeared correlator C{t) with interpo¬ 
lating field J'dd on the lattice of size 32® x 64. The fit range 
is 5 — 20 in lattice units and the values at the chiral limit are 
obtained by linear extrapolation. It has been verified that the 
values for the masses do not change outside errorbars if the 
lower fit range is increased from 5 to 7-10. 


Negative Parity 

Positive Parity 

Kl 

am 

am'^ 

0.153 

1.226(17) 

1.366 (35) 

0.1550 

0.999(21) 

1.191(45) 

0.1554 

0.946(23) 

1.167(48) 

0.1558 

0.891(25) 

1.149(54) 

0.1562 

0.834(29) 

1.135(63) 

f^c 

0.701(29) 

1.036(51) 


mass values leads to the determination of the mass in 
physical units. We find in the negative and positive par¬ 
ity channels the values 

m- = 1.605 ± 0.058 GeV m+ = 2.36 ± 0.13 GeV (22) 
respectively. As we have already pointed out the ratio 
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FIG. 18: The mass squared of the kaon (crosses), nucleon 
(open triangles), pentaquark in the negative parity channel 
(filled circles) and in the positive parity channel(asterisks) 
plotted versus the pion mass squared. We use the nucleon 
mass at the chiral limit to set the scale obtaining a~^ = 
2.2 GeV. The dashed lines show the linear extrapolations to 
the chiral limit. The dotted line shows the KN threshold at 
the chiral limit. 
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FIG. 19: The mass of the lowest eigenstate of the pentaquark 
system in the negative parity channel as a function of the pion 
mass squared. Filled squares show the results of this work, 
crosses are data from Ref. M , open triangles from Ref. .HI, 
open circles from Ref. open squares from Ref. 113 and 
filled circles from Ref. flil . In each case we give the spatial 
length of the lattice. 
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between the mass of the 4> meson evaluated at Kg = 0.155 
and the mass of the nucleon at the chiral limit is 1.002 ± 
0.025, within 10% of the value of 1.087 obtained using 
the physical (j) meson and nucleon masses. Similarly the 
ratio rriK/mN = 0.493(17) at the chiral limit is very close 
to the physical value of 0.526. 

Finally we compare our analysis to other lattice studies 
of the pentaquark system. All lattice studies so far are 
carried out in quenched QCD. The first pioneering stud¬ 
ies with Wilson fermions were carried out by Csikor et 
al. ^3 using a variant of the KN interpolating field, 
where the color index of the s is coupled to a light quark 
in the nucleon instead of in the kaon and S. Sasaki M 
using the diquark-diquark interpolating field. Both found 
evidence for a pentaquark state. Motivated by these 
studies we looked at the pentaquark density-density cor¬ 
relator in order to learn about the internal structure of 
the 0+ [l^ . Chiu and Hsied ^ using overlap fermions 
performed a variational analysis using interpolating helds 
Jkn, Jkn ^ud Jdd- Their conclusion was that there 
is a resonance in the positive parity channel, which in 
the chiral limit has a mass close to that of the 0+(154O). 
The next lattice group to have results on the pentaquark 
system was Mathur et al. 11^ . They used interpolat¬ 
ing helds Jkn and J'j^^ [l^l- By studying the scaling 
of spectral weights on two volume of size 12^ x 28 and 
163 X 28 using a sequential Bayesian method they con¬ 
cluded that the states that they observed were KN scat¬ 
tering states. In the work of Ishii et al. in addition 
to periodic boundary conditions, antiperiodic boundary 
conditions were used in the spatial direction for the light 


quarks. With antiperiodic boundary conditions the low¬ 
est momentum allowed for each quark is tt/L. Therefore 
in the negative parity channel switching from periodic to 
antiperiodic b.c. should increase the energy of the lowest 
KN scattering state. They indeed observed the expected 
shift in energy concluding that the lowest state is a KN 
scattering state. Our analysis is closest to the analysis 
carried out in ref. M- Diagonalizing the 2x2 correla¬ 
tion matrix constructed using Jkn and J'kn^ using 
a large number of conhgurations they were able to accu¬ 
rately determine the two lowest eigenvalues in the neg¬ 
ative parity channel and check scaling of their spectral 
weights. Based on this scaling they concluded that the 
lowest state is the s-wave KN scattering state whereas 
the second is a resonance state that they identified as 
the O'*". In the positive parity channel they obtained 
only one state for which they concluded that it is not 
the p-wave KN scattering state. Finally after the com¬ 
pletion of our work a study of the binding energy in the 
pentaquark system as a function of the light quark mass 
was carried out in Ref. [l^. The fact that they saw no 
increase in the binding as the light quark mass decreases 
led them to conclude that there is no pentaquark reso¬ 
nance. Following Ref. |l 9|| we compare our results given 
in Table for the negative parity channel to the results 
from other lattice groups in Fig. There are small 
variations in the measured data among different groups 
which, at least partly, can be attributed to small differ¬ 
ences in the strange quark mass. Errors due to the choice 
of the lattice scale a are also not included. 
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IV. CONCLUSIONS 

We have presented a computation of the mass of the 
pentaquark system using interpolating fields, which are 
motivated by the diquark-diquark and diquark-triquark 
description proposed for the internal structure of the 0“''. 
Despite the difference in the structure of these interpo¬ 
lating fields the values obtained for the lowest mass of the 
pentaquark system using either interpolating field are in 
agreement with each other. They are also in agreement 
with the lowest energy eigenvalue determined from the 
analysis of the correlation matrices constructed using ei¬ 
ther local or smeared KN and diquark-diquark interpo¬ 
lating fields as a basis. We check our lattice techniques by 
studying the two-pion system in the 1=2 channel where 
no low lying resonances are present. By analyzing the 
correlation matrix constructed using products of pion 
and rho-meson interpolating fields we accurately deter¬ 
mined that the two lowest states are the s-wave two-pion 
and two-rho scattering states with a mass gap which is 
about 320 MeV at a quark mass that corresponds to pion 
mass of about 830 MeV. It is explicitly verified that these 
are s-wave states by projecting to zero relative momen¬ 
tum on our smallest lattice. We note that the energies 
obtained by diagonalizing the correlation matrix in the 
center of mass frame are a mixture of scattering states 
of different relative momentum even when the two low¬ 
est states are scattering states and only at large time 
separations they yield the correct ground state. This is 
a relevant result since in many lattices studies it is as¬ 
sumed that a similar diagonalization resolves the scatter¬ 
ing states with different relative momentum. By studying 
the spectral weights for the two lowest eigenstates in the 
two-pion system on our three lattice volumes we show 
that the correct scaling with the spatial volume sets in 
at large times requiring accurate determination of local 
correlators. The spectral weights in the pentaquark sys¬ 
tem show no volume dependence in the time range where 
in the two pion system we saw a clear scaling of spectral 
weights. Therefore based on our spectral weights results 
we can not exclude a resonance pentaquark state. In 
the negative parity channel our correlation matrix anal¬ 
ysis gives accurately one state close to the KN threshold 
with a first excited state that can be determined only 


at very short time separations. The mass of the low¬ 
est state is consistent with the mass extracted from just 
Jdd and Jkn- We would like to stress that the fact that 
within the time extent of our smaller lattices the spectral 
weights do not scale with the volume only leads to the 
conclusion that we do not have a single scattering state. 
Given the fact that within our statistics and time ranges 
we are unable to accurately resolve a lower KN- scatter¬ 
ing state and a higher close-by single particle state does 
not permit us to draw any definite conclusion regarding 
the existence of the 0'*". However we have shown in this 
work that one would need very accurate data, better in¬ 
terpolating fields and lattices with large time extension 
in order to reliably resolve the low lying states and per¬ 
form a spectral weights analysis to exclude or show the 
existence of the 0'*'. Given that the lowest eigenstate 
from the correlation matrix analysis yields an energy in 
agreement with that extracted with J^dd means that we 
can use the local-smeared diquark-diquark interpolating 
field to compute the mass of this state. This is carried out 
on a set of five values of light quark masses on our large 
lattice. A linear extrapolation to the chiral limit leads to 
a value for the mass which is about 10% above the KN 
threshold. In the positive parity channel the two low¬ 
est eigenstates are separated by an energy gap of about 
100 MeV at the two heaviest pion masses. However the 
mass that we find in the chiral limit, determined in the 
same way as in the negative parity channel is about 920 
MeV above the KN threshold and therefore too high to 
be identified with the 0’'"(154O). In summary, we have 
shown that within this analysis in quenched lattice QCD 
and using Wilson fermions, we cannot exclude a pen¬ 
taquark state, which in the negative parity channel has a 
mass about 750 MeV lower than in the positive channel. 
However, with the lattice sizes used and within our statis¬ 
tics, the existence of the 0-1- has also not been established 
from this study. In order to reach a definite conclusion 
regarding its existence one would require a more detail 
and accurate computation involving a larger and better 
basis of interpolating fields, lattices with larger temporal 
extent and more statistics. 
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